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PACS 71.38.-k - Polarons and electron-phonon interactions 

PACS 72. lO.Di - Scattering by phonons, magnons, and other nonlocalized excitations 
PACS 63 . 20 . kd - Phonon-electron interactions 

Abstract. - We investigate the existence of end points in the dispersion of Holstein polarons in 
various dimensions, using the Momentum Average approximation which has proved to be very 
accurate for this model. An end point separates momenta for which the lowest-energy state is 
a discrete level, i.e., an infinitely-lived polaron, from those where the lowest-energy feature is a 
continuum in which the "polaron" is signalled by a resonance with a finite lifetime. While such 
end points are known to not appear in ID, we show here that they are generic in 3D if the particle- 
boson coupling is not too strong. The 2D case is "critical": a pure 2D Holstein model has no end 
points, like the ID case. However, any amount of interlayer hopping leads to 3D-like behavior. 
As a result, such end points are expected to appear in the spectra of layered, quasi-2D systems 
described by Holstein models. Generalizations to other models are also briefiy discussed. 
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Introduction. — Understanding the properties of a 
novel material is inherently linked to understanding the 
physics of its underlying quasi-particles. The formation 
of such quasi-particles is a very old problem that comes 
up again and again in condensed matter physics. In par- 
ticular, the formation of a quasi-particle qp due to the 
coupling of the particle to bosonic degrees of freedom is 
relevant in many systems of current interest, such as man- 
ganites, cuprates, most other oxides, various organic mate- 
rials, etc [ ~ ]. This coupling can be to phonons, magnons, 
orbitons, or any combination thereof, and the resulting 
composite object is known as a polaron. 

Such polaron problems continue to attract considerable 
analytical and numerical interest, particularly for simple 
lattice Hamiltonians like the Holstein model [ ] vifhich can 
be studied by many means [ ]. The ground-state is a dis- 
crete eigenstate, i.e., it describes a true quasiparticle with 
an infinite lifetime. At higher energies, however, the eigen- 
states may describe incoherent scattering of the particle on 
bosonic excitations, and have finite lifetimes. 

One important question which has not been answered 
for lattice models of polarons, and which we settle for 
the Holstein polaron in this Letter, is whether the low- 
est eigenstate for a given momentum is an infinitely-lived, 
discrete state throughout the Brillouin zone, or not. That 
such a state may not exist for all momenta has long been 
suspected, since the dispersion of the Frohlich polaron ex- 
hibits a so-called end point \ ]. While this is a three- 



dimensional, continuous model, it is interesting to see if 
similar phenomenology is possible for lattice models, and 
under what conditions. 

As we show below, the answer depends strongly on both 
the dimensionality of the system, and the value of the ef- 
fective particle-boson coupling. At first sight, the former 
seems surprising, since many authors have demonstrated 
that dimensionality plays a very small role in determining 
the qualitative features of the ground state properties of 
the Holstein polaron ['^-Id]. However, the bare electron 
and polaron densities of states vary significantly with di- 
mension [11], and as we show below, this plays an essential 
role in deciding whether a discrete level is pushed below 
the polaron-|-one-phonon continuum, or not. This con- 
tinuum appears at f2 (the energy of the boson, modelled 
as an Einstein mode) above the ground-state energy, and 
describes excited states with a boson far from the polaron. 

That the effective coupling is also very important is, 
on the other hand, not surprising. At very large cou- 
plings, the polaron effective mass is exponentially large, 
and therefore the polaron dispersion is very flat. As a 
result, one expects that the entire polaron band fits be- 
low the continuum, i.e., a true quasiparticle state exists 
at low energies in the entire Brillouin zone. At weaker 
couplings, however, the polaron dispersion is considerable 
and its band may overlap with the continuum at larger mo- 
menta. For the most studied case - the ID chain - it is well 
known that irrespective of the coupling, (at least) one dis- 
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Crete state is pushed below the continuum everywhere in 
the Brillouin zone, although it flattens out just underneath 
the continuum at larger momenta. While it seems to be 
assumed in the literature that this is also the case in higher 
dimensions, we show here that this is not true in 3D. In- 
deed, in 3D at weak couplings, an infinitely lived polaron 
exists only in a finite region near the F point. The 2D case 
is marginal: while a low-energy discrete state always exists 
in a purely 2D model, any amount of anisotropic hopping 
in the 3rd direction renders the system effectively 3D-like, 
and the transition from a true quasiparticle to a resonance 
becomes possible if the coupling is not too strong. In other 
words, unlike ID systems, layered quasi-2D systems are 
not guaranteed to have a long-lived quasiparticle as the 
low-energy eigenstate throughout the Brillouin zone. 

Formalism. — Our results are based on the Momen- 
tum Average (MA) approximation, initially developed to 
study the properties of the Holstein model [J n] and then 
extended to more complex polaron models [12]. The ap- 
proximation allows one to sum all diagrams appearing 
in the expansion of the self-energy, albeit with exponen- 
tially small contributions discarded so that this sum can 
actually be performed analytically. The MA approxima- 
tion has been shown to be accurate over all of parameter 
space (except the extremely adiabatic limit), in all dimen- 
sions, and for all energies. Furthermore, the approxima- 
tion satisfies the first six spectral weight sum rules exactly, 
and all higher order sum rules to a good degree of accu- 
racy. It can also be systematically improved by including 
additional states, allowing for the correct reproduction 
of the polaron-|-one-phonon continuum. These improve- 
ments also lead to increased accuracy throughout param- 
eter space, and exact agreement with additional spectral 
weight sum rules. Recently, the MA approximation has 
also been generalized to treat broken translational sym- 
metry, due to disorder or surfaces [1.3], and also to ob- 
tain two-particle Green's functions needed for a calcula- 
tion of the optical conductivity [14]. In conclusion, the 
MA approximation is a powerful and well-understood tool 
that allows us to understand very accurately the Holstein 
polaron physics, throughout the parameter space, for all 
energies and momenta. Also, because it is an analytical 
approximation, we can easily probe large regions of pa- 
rameters space that are difficult for numerical techniques 
to investigate, especially in higher dimensions. 

The Holstein model that we analyze here describes the 
simplest possible electron-phonon coupling on a lattice, 
and is given by [']: 



simple hypercubic lattice of lattice constant a: 
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The coupling is to a branch of Einstein optical phonons 
of frequency ft, where and are the phonon cre- 
ation and annihilation operators. The last term de- 
scribes a momentum-independent on-site linear coupling 
(7^jc|cj(5| + 6j), written in k-space. All sums over mo- 
menta are over the first Brillouin zone. For all results 
shown here, the total number of sites N is taken to be 
infinity. We also set h = 1 and a = 1 throughout. 

The MA approximation has been applied to the Holstein 
model previously, and its meaning and accuracy is very 
well understood [ ]. Although the MA*^°^ level of approx- 
imation has been shown to be quite accurate at describing 
ground state properties, it fails to reproduce the correct 
location for the polaron-|-one-phonon continuum. At the 
MA^^) level and beyond, this continuum is accurately re- 
produced. Although the discussion and conclusions below 
do not qualitatively depend on the precise location of this 
continuum, we will employ the more accurate ap- 
proximation throughout this work so that our results are 
quantitatively, not just qualitatively, accurate. 

Application of the MA approximation to the Holstein 
Hamiltonian allows us to calculate the single-particle (re- 
tarded) Green's function, defined as G'(k,w) = (0|c|^(a; — 
V, + ii])^^ cl^\0) [15]. The solution for this Green's function 
can be written in the standard form [lU]: 



G(k,w) = (2) 

w - £k - (k, w) + irj 

with the MA'-"'^^ self-energy given by: 

y /, ,^ g'go('^) 

1 - ggQ(uj) [A2(uj) - Ai(w - i2)J 



where 



50 (w) = j^^^GoiK 



(4) 



is the momentum average over the Brillouin zone of the 
free propagator G'o(k,CL') = {uj — e\^ + i'r])~^ . The functions 
An{uj) are infinite continued fractions, defined as: 



1 - 55o(w - nn)An+i{u}) 

nggo{uj - nfl) 



(5) 



1 



{n + l)g^go{uj - nn)go{u} - {n + l)f2) 



1 



The first term is the kinetic energy of the free electron, 
with being electron creation operators for a state with 
momentum k (the spin is irrelevant in this model, and we 
ignore it in the following). For the free electron disper- 
sion, we use nearest-neighbor hopping on a d-dimensional 



We have also used the short-hand notation ui = to — ft — 
gAi{uj—fl). As discussed at length elsewhere [10], this MA 
self-energy is k-independent simply because the Holstein 
model is so featureless. Indeed, momentum dependence 
for E(k, w) is obtained for the Holstein model for MA*^") 



p-2 



The end points in the dispersion of Holstein polarons 



with n > 2 (although it is very weak). For other mod- 
els with momentum-dependent couplings, a k-dependent 
self-energy is found even at the MAf^"^ level [12]. Further 
details on the derivation and meaning of the MA approx- 
imation can be found in Ref. [ ] . 

With an explicit expression for the Green's function of 
the Holstein polaron, we calculate the spectral weight 

^(k,w) = --ImG'(k,w), (6) 

TT 

from which we can extract the polaron dispersion, av- 
erage number of phonons, effective mass, etc. [10]. In 
our calculations we employ a small but finite value for 
J]. This moves the poles of the Green's function off of the 
real axis and changes the (5-peaks of the spectral weight 
A(k,uj) = J2a l(Q!|Ck|0)|^<5(i^ — Ea) into Lorentzians. In 
practice, it is necessary to choose 77 small enough to allow 
detection of the Lorentzian peaks in regimes where the qp 
weight is extremely small. It is also convenient to define 
the dimensionless coupling constant A = g^/{2dVlt). We 
now turn our attention to the conditions that lead to the 
formation of an infinitely-lived polaron. 

Results. — Separating the self-energy in terms of its 
real and imaginary parts, — + we can 

write the spectral function in the following form: 

^ 1 ^-s"M .7. 

^ ''^^ ^[cj-ek-S'(u;)]2 + [,7-S"H]2- ^' > 

We are interested in the dispersion of the polaron, E\^. 
This can be found by tracking the lowest pole in the spec- 
tral weight [10], its energy being the smallest value of lo 
which satisfies 

w-Ek-S'H. (8) 

For this to also be a discrete eigenstate (as opposed to 
the lower edge of a continuum) , the imaginary part of the 
self-energy has to vanish: 

S"(^)l.=^,,-0. (9) 

The latter requirement is essential, as eq. (7) does not 
have singular behavior when the imaginary part of the 
self-energy is finite. If eq. (9) is satisfied, we call the 
solution a true polaron qp, meaning that it has an infinite 
lifetime. If eq. (9) is not satisfied, then the lowest-energy 
feature is not a discrete state but the polaron-|-one-phonon 
continuum. In this case, we find the energy where this 
continuum has a maximum in the density of states (DOS), 
and call it a polaron resonance. This resonance has a finite 
lifetime proportional to 1/E"(a;), where ui is its energy. 

As already mentioned, the onset of the continuum is 
at Eq + ri, where Eq is the ground state energy of the 
polaron for k = [10], therefore for a; < £^0 + ^ it is 
guaranteed that 5]"(a;) = 0. To summarize, we thus look 
for the region(s) of parameter space where i^k, the lowest- 
energy solution of eq. (8), satisfies E\^ < Eq -\- Q,. If this 
condition is satisfied, we have a true polaron bound state 
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Fig. 1: (Color online) (a) The real part of the ID self-energy 
E(cj) versus u)/t, and (b) the polaron dispersion Ek, for effec- 
tive el-ph coupling strengths A = 0.1,0.3 and 0.6. In (a) we 
also plot uj — Ek, for fc = (solid blue line) and k = tt (dashed 
blue line). The intersection of these curves with ReS(a;) give 
the polaron energy eigenvalues, shown by the respective cir- 
cles in (a) and (b). The shaded areas indicate the onset of the 
continuum, i.e., where lmE(ci;) 7^ 0. 

at that momentum. Otherwise, we have a polaron-like 
resonance with a finite lifetime lying inside the continuum, 
as described above. 

Obviously, it is the self-energy that controls which of the 
two possible cases is encountered. To understand how this 
occurs, we first quickly analyze the well-known ID case, 
where - as already mentioned - there is always a true po- 
laron qp state below the continuum. The reason for this 
is made obvious by fig. 1(a), where we show a geometric 
solution of eq. (8) for three different values of the effective 
coupling A. The real part of the self-energy, S'(w), goes 

2 

asymptotically like ^ as a; — > — oo, so it starts as a nega- 
tive function which, in ID, diverges and changes sign at a 
value marking the onset of the polaron-|-one-phonon con- 
tinuum (above this energy, in the shaded regions, S"(a;) 
becomes finite). The left-hand side of eq. (8) is a straight 
line - the limiting A; = and k = n are drawn explicitly - 
and a true polaron solution exists if the two intersect be- 
low the onset of the continuum. The divergence of S'(a;) 
ensures that in ID there is always such a solution. Its 
energy i?k is shown in fig. 1(b), where we also plotted the 
onset of the continuum, Eq + il, for the three cases. 

As expected, at very weak couplings and sufficiently 
large k, i?k fiattens out just below the continuum. Ge- 
ometrically, we see that this is due to the fact that for 
a large region of the Brillouin zone, the intersection be- 
tween the two curves occurs very close to the divergence 
in S'(w), so the value of the solution ui varies little and is 
asymptotically close to Eq + il. The value of k where this 
"flattening" begins is also easy to estimate for A — 0: in 
this limit, to zero order in perturbation theory we must 
have -Ek ~ £k, therefore the "flattening" must happen for 
momenta k such that ei^ > Eq + ^l. Here, because of the 
divergence in the free-electron DOS at the band-edge, ar- 
bitrarily weak coupling to the phonons suffices to repel a 
bound state below the continuum. As the coupling A in- 
creases, the self-energy curve shifts to lower energies, but 
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Fig. 2: (Color online) (a) The real part of go{ijj) and (b) the 
bare electron DOS (— l/7r)Imgo('^), shown for d = 1,2,3. A 
broadening factor of r]/t — 0.005 has been used. 

the divergence is always there. 

Clearly, the next step is to understand why the self- 
energy has this shape and in particular the reason for the 
divergence in its real part. If such a divergence is always 
guaranteed to mark the onset of the continuum, then a 
discrete polaron state will always be found. 

As already discussed, we are actually primarily inter- 
ested in the limit of weak el-ph coupling, since at large 
el-ph coupling the polaron band is very flat and guaran- 
teed to obey i?k < Eq + il everywhere in the Brillouin 
zone. For weak coupling g — > 0, eq. (3) gives a self-energy 



i;(a;) « g^goi^ - fi) + 



N 



after using eq. (4). We recognize this as being the lowest 
order diagram contributing to the self-energy (the Born 
approximation), in agreement with perturbation theory. 
The correction w — — >■ cj in the argument of go('^) will 
further shift the self-energy towards lower energies, and is 
due to corrections from higher order diagrams. However, 
this shift cannot be responsible for the appearance of a 
singularity in the real part of the self-energy. 

It follows that for the Holstein model (we discuss other 
models below) the real part of go{oj) controls the energy 
dependence of 5]'(cj) for small A, and therefore whether it 
has a singularity or not. Both the real and imaginary parts 
of goiijj) are shown for various dimensions in fig. 2. Ac- 
tually, the right panel shows (— l/7r)Imgo('^), which from 
eq. (4) is seen to be the bare electron total DOS. 

Because go('^) is the momentum average of a retarded 
propagator, its real and imaginary parts are related 
through Kramers-Kronig relations. Given that the imagi- 
nary part (the DOS) is finite only inside the free-electron 
bandwidth < 2dt, it follows that it is the DOS at the 
band-edge that controls whether the real part has a sin- 
gularity at the band-edge energy: if the band-edge DOS 
diverges, or at least has a discontinuous jump, then the 
real part has a singularity at the band-edge; otherwise it 
is finite everywhere. 

This explains why in ID we always find a singularity in 
the real part of ga(u!), and therefore in S'(w): it is well- 



Fig. 3: (Color online) (a) The real part of the 3D self-energy 
for an intermediate el-ph coupling A — 0.7. The shaded region 
shows where Im E(tj) 7^ 0. The straight lines show cj— Sk curves 
for ki = 7r(0.2, 0.2, 0), ka = 7r(0.3, 0.3, 0), kg = 7r(0.43, 0.43, 0). 
An intersection of oj — and Re corresponds to a solution 
of eq. (8); (b) Contours of the end points for various el-ph 
coupling strengths. The A = 0.7 contour is shaded, and the 
arrow marks the end point kc = ka found in (a). 



known that the free electron DOS for nearest-neighbor 
hopping in ID has van Hove singularities at the band 
edges. In 2D, the DOS has discontinuities at the band 
edges. As a result, based on the discussion above, we ex- 
pect that in 2D a discrete polaron state must also exist in 
the entire Brillouin zone, just like in ID. However, this is a 
very weak logarithmic singularity (as it is due to a discon- 
tinuity, not a singularity in the DOS) and one may expect 
that it can be easily removed. This is indeed the 
we show below where we analyze quasi-2D systems. 

However, first we focus on the 3D case, where the DOS 
at the band-edge is continuous, leading to a real part of 
5o(w), and therefore low-couphng S(a;) = g^goC^), that is 
finite everywhere. The geometric solution of eq. (8) in this 
case is illustrated in the left panel of fig. 3, for an interme- 
diate effective coupling A = 0.7. The low-energy part of 
T,'{oj) still resembles the real part of go{uj), although due to 
higher order diagrams, the shape is somewhat distorted. 
The differences are more considerable at higher energies 
(this issue is revisited below). The onset of the contin- 
uum at Eq + is marked by a discontinuity in the slope 
of S'(a;) (shaded area). Because T,'{uj) does not diverge 
below this energy, it is now apparent that eq. (8) only has 
solutions with E]^ < i?o + up to a critical value kc which 
marks the "end point" in the polaron dispersion ["]. For 
A = 0.7 and along the (1,1,0) direction, the critical value 
is at kc = 7r(0.43, 0.43, 0), as shown in fig. 3(a). 

Repeating this along various directions allows us to 
identify the surface described by the "end points" kc in 
the Brillouin zone, which marks the separation between 
the region near the F point where an infinitely-lived quasi- 
particle exists, and the region near the edges of the Bril- 
louin zone where only a resonance with a finite lifetime 
appears. The end points contours in a quadrant of the 
kz — plane are shown in fig. 3(b), for various values of 
A. The arrow marks the particular kc found in fig. 3(a). 

For a given value of the phonon mode, the smallest sur- 
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Fig. 4: (Color online) The (a) real and (b) imaginary parts of 
the 3D self-energy, for increasing el-ph coupling strengths. The 
self-energy becomes divergent for large enough coupling. 

face of end points is for A — > 0. By analogy with a previous 
discussion, in this case we can estimate it to be given by 
£kc ~ Eo + ^1 since for A — )■ 0, i?k ~ ^k- It follows that 
for any > 0, there is a vicinity of the F point where a 
true polaron exists. In general, however, there are regions 
in the Brillouin zone that do not satisfy this condition 
for small A. Indeed, since the maximum value of Ck is 6t 
while Eq = — 6i, we estimate that for any f2 < 12t there 
exist such regions of momenta without an infinitely lived 
polaron state, < 12t is likely to hold in most materi- 
als, therefore such a case is rather typical. If J7 is not too 
large, then in the limit A — the surface of end points 
is spherical, since at low energies the free-electron disper- 
sion can be approximated by a parabola. As either $1 or 
A increase, the shape will expand and distort, as seen in 
fig. 3(b) for increasing A. For any value of < 12t, there 
is a critical value Ac > above which an infinitely-lived 
polaron solution appears in the entire Brillouin zone. For 

= 0.5t, fig. 3(b) reveals that this Ac is just above 0.86. 

As already discussed, the fact that a true polaron solu- 
tion exists everywhere in the Brillouin zone at strong(er) 
coupling is not surprising, because the significant increase 
in the polaron effective mass guarantees that eventually 
the whole polaron band will fit below the continuum. In 
terms of the geometrical solution, what happens is re- 
vealed in fig. 4 which shows that for large enough A, the 
real part of the self-energy does eventually gain a singu- 
larity at low-energies. This divergence is a higher order 
effect, not coming from a divergence/discontinuity in the 
DOS, but due to the structure of the self-energy itself. 

Eq. (3) reveals that at large enough el-ph coupling g, 
the denominator itself can vanish giving rise to a different 
type of divergence in the self-energy. The appearance of 
this divergence explains the phenomenology related to Ac 
(for 51 = 0.5i, this divergence is first observed at Ac — 
0.87). Physically, it is linked to the formation of the second 
bound state [10, 16], which is known to form below the 
polaron-f one-phonon continuum as the crossover towards 
the small polaron regime is approached for A ~ 1. 

For completeness, we mention that as the coupling in- 
creases further, more and more divergences appear in 
E'(a;) due to the structure of the continuous fractions, see 
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Fig. 5: (Color online) (a) The spectral weight A(k, uj) versus 
ui with ky = = 0, and k^/n = 0,0.2,0.4,0.6,0.8,1.0. The 
dashed line denotes the onset of the continuum. With the 
parameters shown, there is no longer a polaron qp for k^/n > 
0.46, but a polaron-like resonance, (b) The lifetime of the 
polaron-like resonance as a function of momentum. The dashed 
lines denote the transition to an infinitely lived polaron state. 

eq. (5). These additional divergences occur at higher and 
higher energies that are well above the onset of the contin- 
uum, thus they are not relevant for the problem we con- 
sider here (they are linked to the crossover towards a Lang- 
Firsov type of spectrum expected in the limit A — >■ oo). 

What happens for A < Ac as we sweep through kc 
is shown in fig. 5(a), where we plot the spectral weight 
A(k,w) vs. u!, for several values of k. The onset of the 
continuum is marked by the dashed vertical line. For mo- 
menta below kc, we see a Lorentzian peak below the con- 
tinuum. Its width is controlled by rj. This is the discrete 
level whose energy is i?k- As k increases, i?k approaches 
the continuum. Unlike in ID and 2D cases, however, where 
the singularity/discontinuity in the DOS ensured that a 
discrete level is always pushed below the continuum, here 
we see that above kc, the low-energy feature is the con- 
tinuum, with a broad maximum that disperses towards 
higher energies while broadening even more as k increases 
- this is the polaron resonance. Its finite lifetime, equal 
to 1/S]"(a;) where uj is its energy, is displayed in fig. 5(b) 
along various cuts in the Brillouin zone. As expected, it 
diverges as kc is approached from above, and it decreases 
very fast towards the edges of the Brillouin zone. This is 
a T = calculation, so this finite lifetime is an intrinsic 
effect. As a technical note, in calculating it we decrease rj 
until changing it any further has no effect on the lifetime. 

Based on all these results, we can now quickly analyze 
the 2D case. For a purely planar 2D hopping model, the 
well-known discontinuity in the DOS at the band-edge 
guarantees that a low-energy singularity always exists in 
I]'(a;), see fig. 2. As a consequence, a true polaron is guar- 
anteed to form for all k at any A, just like in ID. How- 
ever, most layered materials are in fact quasi-2D, because 
there is some small anisotropic hopping between different 
planes. The effect of a small — O.lt on the go{oj) func- 
tion is shown in fig. 6. Any ^ will smooth out the 
band-edge discontinuity in the DOS, and thus remove the 
weak singularity in the real part of go('^) and therefore 
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Fig. 6: (Color online) Anisotropic hopping - transition from 
2D to 3D DOS. 

in S'(a;) at weak coupling. As a result, quasi-2D systems 
behave similarly to 3D systems: for weak-to-intermediate 
coupling, we expect that a surface of end points separate 
the region near F where a true polaron quasiparticle forms, 
from a region near the Brillouin zone edges where only a 
finite lifetime resonance appears. The detailed shape of 
this surface and the critical couplings Ac above which a 
true polaron exists everywhere will, of course, depend on 
the particular parameters of the problem. However, the 
reasons we uncovered here to explain this phenomenology 
are robust and therefore such surfaces of end points should 
be quite generic in quasi-2D and 3D materials with not too 
large el-ph coupling. 

The existence of such surfaces of end points, especially 
in layered systems, should be directly visible using Angle 
Resolved Photo-Emission Spectroscopy [ ]. This tech- 
nique measures the spectral weight and should be able to 
identify significant changes in its shape, like those shown 
in fig. 5(a). Very different temperature dependence of fea- 
tures in different regions of the Brillouin zone would also 
be consistent with such phenomenology, because the true 
qp peaks should only exhibit thermal broadening, while 
the finite-lifetime resonances have an intrinsic lifetime. Of 
course, many other properties should be strongly influ- 
enced by the existence of such a surface of end points, 
especially if the Fermi surface is not too far from it (this 
statement assumes that interactions between polarons, ei- 
ther due to direct interactions between the particles or 
mediated through exchange of bosons between the two 
polaronic clouds, do not qualitatively change this phe- 
nomenology. Whether this is the case or not is a very 
interesting question which still awaits resolution). 

To summarize, for a Holstein model and weak-to- 
intermediate coupling one can deduce easily whether an 
infinite lifetime polaron qp is expected everywhere in the 
Brillouin zone or not simply from knowing whether the 
DOS of the free electrons is continuous or not at the 
band-edge. To pinpoint the surface of end points, more 
detailed calculations like the ones shown here are needed. 
For more complicate polaron models, with a particle-boson 
coupling that depends on either (or both) the particle and 
boson momentum, and also with a possibly dispersing bo- 



son mode, one can certainly gain intuition about what is 
likely to happen from studying the Born approximation 

S(k, ^) « ^ E l3k,qpGo(k - q,c. - r!q) 
q 

and whether it generically predicts infinitely-lived quasi- 
particles for all k or not (if there are no end points at 
weak coupling, they are very unlikely to appear at stronger 
couplings for reasons already discussed). If either gt.q or 
flq have strong momentum dependence, then S(k, w) may 
acquire a sufficiently strong k-dependence to modify the 
conditions found above for the Holstein model. In particu- 
lar, it is not apriorily obvious that a singularity in S'(k, uj) 
is absolutely needed to guarantee a discrete low-energy so- 
lution everywhere in the Brillouin zone, if its momentum 
is strong enough to shift the curve around significantly. 

^ ^ ^ 
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